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I. INTRODUCTION 



The concept of quantum non-demolition (QND) measurement was originally introduced 
in the context of the detection of the gravitational waves [|I| which produce very weak 
changes in the measuring instrument. In general, the quantum mechanical unpredictability 
disturbs the system being measured. This presents particular challenges in the context 
of ultra-sensitive measurements where the disturbances may mask the real effects to be 
measured. A sequence of completely predictable measurements on a quantum system, in the 
sense that the measurements do not disturb the system ("back-action evading"), is, however, 
possible for certain detector observables A{t), called QND observables, that satisfy the 
relation [^(t), A{t')] = 0, where t and t' denote different times. In addition to its relevance in 
ultra-sensitive measurements, a QND scheme provides a way to prepare quantum mechanical 
states which may otherwise be difficult to create, such as Fock states with a specific number 
of particles. In the context of research in Bose-Einstein condensation (BEC), such quantum 
state preparation procedure would be particularly valuable, for instance, in engineering 
entangled state or Schrodinger's cat states. Indeed, this is perhaps the most important 
motivation in searching for QND schemes with atoms. It is noted that a state preparation 
with BEC has recently been performed in the form of squeezed state creation in an optical 
lattice 10. A successful QND scheme would create an important quantum state in the form 
of the Fock states of known number of atoms. 

One of the original proposals for a quantum optical QND scheme was that involving the 
Kerr medium which changes its refractive index as a function of the number of photons 
in the "signal" pump laser. We wish to apply the same strategy with the atoms where 
the nonlinearity is intrinsic i.e. due to atomic interactions. The advent of experimental 
methods for producing BEC should, in principle, enable us to make progress in the matter- 
wave analogue of the optical QND experiments. In particular, the presence of the long 
range order in a BEC implies that the condensate phase can be accurately measured. The 
schemes to be proposed are based on experimental technologies available immediately: (1) 
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the Bragg scattering of the condensates [Q, and (2) creation of multicomponent, trapped 
Bose condensates in different hyperfine states P,^. 

II. QND USING BRAGG SCATTERING 

Consider a QND scheme in which two hght fields are used to produce equal parts of a 
trapped BEC which is Bragg scattered into distinct momenta hkp and —hkp in the center- 
of-mass (CM) frame |Q. The momentum hkp may be any multiple of unit momentum kick 
hkn, hkp = nhkn, where kji is the wave number of the Raman laser. The relevant parts of 
the field operator at the end of the Bragg pulse may be expanded as a superposition of two 
counterpropagating modes of momenta ±hkp in this CM frame as 



where V is the volume of confinement, and dp, d_p denote the annihilation operators of the 
two counterpropagating modes. Evolution of this Bragg scattered state under the standard 
many-body boson Hamiltonian provides interatomic collisions which are necessary for the 
QND measurement. 

The effective Hamiltonian for this system is then given by 

h^k^ 

^Eff = lTir7'(^p^p + ^-p^-p) + U+pdld^dpap + U-pd'ipdLpd^pd^p + U^pdldpdLpd^p. (2) 

The Hamiltonian is similar to the well-known quantum optical version due to a Kerr medium, 
except for the self-interaction terms i.e. the 3rd and 4th terms. We have deliberately 
distinguished the self-collision and the interaction between the +p and —p modes by using 
U^p, U-p, and U±p. We shall discuss this in more detail below. The first four terms represent 
two self-interacting condensates moving towards each other with momenta —hkp and hkp 
respectively, while the last term provides coUisional coupling between the Boson number 
operators dj^dp and dLpd^p, mediated by the interatomic collisions. We shall refer to this last 
term as the interaction Hamiltonian, Tij. We shall assume that any quasiparticle excitations 
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created by the collision are minimal by assuming that we are in a slow moving regime 
{kp — * 0) and that the energies deposited by the collisions (as determined by the kinetic 
energy corresponding to the relative momenta, 2hkp) do not match the energy of collective 
excitations of a confined condensate. 

The collisional interactions are more accurately described by a representation of the two 
body T-matrix; U+p, U-p and U±p should therefore be calculated accordingly. In particular, 
in modelling a specific experiment one has to treat the expansion of the condensate as 
well as the collisional loss of atoms "knocked out" by other passing atoms in order to 
assess the quality of QND schemes. A "spherical shell" of atoms is generated when two 
distinct momentum wave packets collide, due to elastic scattering into the many momentum- 
conserving channels; a significant portion of the atoms may consequently be lost from the 
characteristic mixing process. It should be possible to perform the whole measurement 
within a trap so that elastic losses are minimised. Recent calculations indicate that the 
fraction of atoms retained in the original wave packets approaches 90% if the asymmetry 
of the trapping potentials is increased to an aspect ratio {uz/uj^ = Uz/uy) of order 0.35 
0]; we assume quasi 1-dimensional systems such as the cigar-shaped condensates in this 
work. Using the result of Ref. [0, we write U+p = U_p = where Uq = and 

f/±p = 2f/o(l + 27(M)- 

The imaginary part of U±p generates elastic losses proportional to the relative momentum 

of the wave packets for higher momenta. In this limit 7(A;) = 2kpa where a is the interatomic 

scattering length; this implies that a slower probe moving through the target is needed to 

minimise the elastic collisional losses. The elastic collisional loss is further reduced by 

superfluidity effects, and energies in the phonon regime {kp — > 0) could be used, as in this 

limit one has 7(fc) oc fc^. In addition, with an initial estimate of the number, the losses can 

be calculated with high accuracy and we shall therefore assume the calculation is carried out 

in an iterative process where the number of atoms measured is fed back into the calculation. 

We note that the interaction strength also changes in reality as a function of position so 

that the coupling in general depends on the shape of the overlapping beams as well as the 
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composition of the counterpropagating wavepackets. It is reasonable to consider ideal cases 
of uniform and real valued effective interaction in this work, with U±p = U^p by virtue of 
kp ^ and the factors discussed above. 

In the absence of coUisional loss, aj^dp and both constants of the motion for 

this system, and hence any function of dl-pd±p may be chosen to be the QND observable. 
We choose to label Os = the QND "signal" variable, and identify the probe variable 

to be the standard Hermitian phase quadrature operator. Op = i {d^p — a^_^ where d-p 
is such that d-p = \Jd^_pd-p+ lexp(i(^_p), with (p_p denoting the phase operator for d_p. 
Further necessary conditions such as [Op^lhii] ^ and [Os^iii] = are clearly satisfied. Hi 
is independent of the phase operator for the signal; this ensures the motion of Os do not 
become unpredictable due to the uncertainty imposed on the conjugate variable of Os- 

The Heisenberg equations of motion for the Boson annihilation operators dp and d-p are 
given by the coupled nonlinear operator equations 

i—f' = ^pO'p + Uodpdpdp + 2UQa^_pd_pdp (3) 
^'~dt' ^ ^p'^-P ^ UooLpd^pd^p + 2Uod^pdpd-p, (4) 

where Up = and TiUq = Uq/V. In general, exact solutions to nonlinear operator equations 
cannot be found. However, we have dl.pd±pd^p = N±pd^p{t), where N±p = dl.pd±p is a con- 
stant in time which enables us to transform to an interaction picture d_p — > d_pe~**-'^*'~^^°^"*'''* 
to write Eq. as 

.dd- 



i- 



^ = 2Cpa_p, (5) 



dt 

where Cp = Uodj^dp is again a constant of the motion. The solution is simply a_p = 
exp(— ■i2Cpt)d_p(0), where t is the duration of the collisional interaction between the "sig- 
nal" and the "probe" condensates. A measurement of the probe phase therefore provides a 
QND information on the signal atom number for a given atomic species. For ^^Rb samples 
{2Cpt) ~ 10 mrad, with iVp ^ 5 x 10^ and t ~ 1ms. One obtains ^ ^ 2 x IQ-H mrad 
per atom. The increase in sensitivity achievable by increasing t may be compromised with 



the extra collisional losses in the system over time. This can be partially, but not fully, 
eliminated by modelling the loss process. 

As a comparison, we note that a related way of performing a QND measurement would 
involve scattering single atoms off a Bose condensate, as a way to probe the cloud without 
disturbing it. We note that in the quantum mechanical framework, a particle of momentum 
k incident on a potential V{x) is analogous to light wave propagation in a medium with 
some refractive index . In a way reminiscent of a pumped Kerr medium, BEC then forms 
effectively a potential barrier as a function of the number of atoms within the condensate. We 
assume the condensate-induced potential is repulsive, although there are of course attractive 
cases. If one assumes that the energy of the incident atoms is, say, twice the peak energy 
of the condensate, the sensitivity of the phase shift, is calculated to be of the order of 
10 mrad per atom for a typical atomic condensate with = 5 x 10^ 0. A typical QND 
scheme may be visualised essentially as a Mach-Zehnder interferometry where the "phase 
shifter" in one arm is the target (or "signal" ) condensate, and the relative phase shift of the 
probe is measured with respect to a reference condensate split off at the first beam splitter. 
The phase of any given wave function is the sum of its spatial and temporal phases and in 
this particular example the pertinent phase shift which encodes the number is spatial. 

We have also studied the dynamics of two colliding condensates have been simulated 
using a 1-dimensional Gross-Pitaevskii Equation (GPE). We have studied the evolution of 
two condensates within a harmonic trap in which one of the condensates is displaced from 
the center while the other is kept at the origin. This initial condition gives rise to a collisional 
interaction between the two in which the displaced condensate passes through the stationary 
one as a part of its simple harmonic motion. As the two condensates overlap at time t = 
they form an interference pattern. At time t = ^ the condensates again evolve into two 
distinct Gaussian-like profiles, although the original shapes are not restored. The dynamics 
of a reference condensate was calculated by carrying out an identical simulation but without 
the stationary condensate. The phase shift at time t = - was obtained numerically. It 
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is noted that at t = ^ the spatial derivative of the phase profile for both the reference 
and the colliding condensate is zero as the condensate momentarily becomes stationary at 
half period. The simulation was repeated for a number of different values of the nonlinearity 
C = NUq. a linear dependence of the phase shift is observed in accordance with the analysis 
above (Fig. 1). 

A possible experimental scheme would involve the techniques of Mach-Zehnder Bragg 
interferometry using a sequence of local f-vr-^ Bragg pulses [0. A similar experiment has 
also been performed at NIST where the trap was kept on. For a QND experiment a |-pulse 
in the place of a 7r-pulse may be applied on one arm of the interferometer. This would, in 
effect, generate a signal condensate out of the probe in one arm, and the number of atoms 
that has been "pushed out" may be ascertained from the resulting changes in the interference 
pattern. The total time it takes to measure the probe, tprobe should be shorter than the time 
it takes for the condensate to evolve via self-interaction: tprobe m/NAnha ~ 0.1s for ^''Rb. 
From the reported experiment [|10|, the overall temporal and spatial scale of the experiment 
would be of the order r ~ 50yus (much shorted than the self evolution time) and L ~ 1mm. 

III. QND SCHEMES IN MULTICOMPONENT BEC 

A system of condensates in multiple hyperfine states is another candidate for QND mea- 
surement, and allows a potentially more flexible experimental geometry. Two overlapping 
condensates of Rubidium atoms in F = l,mF = —1 and F = 2, mi? = 2 states have already 
been produced by the JILA group which allows coherent transfer of atoms from one state 
to another using appropriate sequence of radio-frequency (RF) and microwave transitions. 
On the other hand, the MIT group created three co-existing condensates in the ground 
state of sodium F = l,mp = —1,0, 1 by the use of an optical dipole trap [^]. In order to 
take into account the condensates in different hyperfine levels, the trapping potential and 
the interatomic potentials need to be written in a more general form: Vtrap = J2i V{^)'Pi 
where Vi = is the projector which projects a state into hyperfine level i, and for the 
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interatomic interaction 5{ri — V2) J2j UjVj where Uj = Anajh'^ /m. The Hamiltonian for the 
multi-component system is then 

H = Ha + Hb + Hint, (6) 

where Ha,b are the usual many-body Hamiltonians but with the field corresponding to atoms 
in magnetic levels A, B. The interaction Hamiltonian may be written in terms of mode 
operators a and h in the two mode (zero temperature) approximation: Hint = Woa^afo^fe. 
The atomic states are assumed to be stable in the sense that they do not change their internal 
states on collisions. The QND coupling term between the modes a) a and 6^6 then gives rise to 
identical mathematical features as above. The Hamiltonian for the triplet ^^Na system can be 
shown to couple atoms in all three states so that terms proportional to al]^a_i(a|ai + Soao), 
+ a\ai), and a\ai[a}_ia_i + ajao) appear in the Hamiltonian. Three possible 
QND signal/probe combinations may be exploited in this case, providing extra fiexibility in 
the experiments. The detection of the probe phase requires separation of the two modes, 
which may be achieved by a Stern-Gerlach spin separation (or an RF output coupling). After 
the separation the phase may be measured e.g. by interfering with a previously prepared 
reference condensate. 

We draw the attention to the fact that there are a couple of possible variations to the 
multicomponent QND scheme. One is the case where condensates in different hyperfine 
states are prepared in separate traps i, i = 1, 2, which are then physically overlapped to 
produce collisional interactions between the condensates. This type of procedure has pre- 



viously been proposed in the context of quantum computation for individual atoms ||Tl| , p!2 
It is a more practical proposition for atoms held in optical traps, whose position can be 
manipulated with relative ease. The Hamiltonian is of the form 

i=i j=i ij=i 

where H^Ai^) Hamiltonian for the component A{B) in the trap with trajectory r^^^\t), 

given by the single particle Hamiltonians plus the nonlinear interactions. The probe and the 
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'H = j:H\ + j:K+j: U,, J dr¥}ir,t)¥^ir,t)¥^ir,t)¥Ar,t), (7) 



signal variables are identified by the different trap/hyperfine state combinations a,f3,A,B. 
Following collision- induced phase shift, the measurement of the phase may be carried out by 
either turning off the trap potential that contains the probe condensate or alternatively, by 
first displacing the trap with the probe. This is schematically illustrated in Fig. 2. Another 
possibility is where the traps of different trajectories contain two-species condensates in 
each. Phase changes in the atomic wave function may be induced such that two atoms in 
different hyperfine states A and B from traps with trajectories r°'(t) and r^(t) undergo a 
transformation ip'Xi'B ~^ ^'''^'^'^ i^li^B where (j)^^ is the collisional phase shift [|l^. For the 
case where the collisional phase shift (p^^ is vr, one may write, using the Schwinger angular 
momentum operators J*, J*, the total Hamiltonian as 

n = QJ: + fi^jf + (j;)' + [J^y + R^pJ^Jl, (8) 

where VL, Ra, Rp, Rap are constants. In this case, both J" and Jf , the particle number 
difference of the two species in traps a, /3, are constants of the motion. The equations of 
motion for the angular momentum raising operator = d\hi in the trap of trajectory i are 
given by 

= n^Jl + 2RaJ':Jl + RpcJ^zJl (9) 

dp 

= npjl + 2RpjPJl + RapJ^Jl. (10) 

The particle number difference J° in trap a may therefore be deduced from measuring J^, 
the phase difference between components and hi in trap (3. J+ may be measured by 
passing the two condensates through an atom beam splitter and measuring the intensity at 



one of the output ports |T3[. All that is required is a way to couple out from the trap a 
superposition of the two components, di + e*'^''='6j, where 0rei is the relative phase between 
the two components. To this end, one may first separate di and hi from the "probe" trap 
into two separate traps (this is possible because they are in different hyperfine states) and 
interfere them normally by opening the traps while "eliminating" any information about 
from which trap an atom came from. This produces a spatial interference pattern which 
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is exactly analogous to the case of having a beam splitter and an atom detector where the 
intensity at one of the output port oscillates in time (rather than space). 

IV. DISCUSSION 

We have discussed several scenarios which enable QND measurement using the experi- 
mentally produced BEC in alkali atoms. The schemes proposed are already within current 
experimental capabilities, and should be achievable in the nearest future. They represent 
a new addition to the currently growing repertoire of atomic quantum optical experiments 
based on BEC. A recent experimental development that requires attention is the achieve- 
ment of BEC in metastable He [Q; its potential sensitivity at the single atom level would be 
suitable for checking QND experiments measuring exact number of atoms. A possible future 
application of the QND scheme is in quantum computation with BEC; it is essential to know 
the exact number of atoms in the condensate in order that the computational operations can 
be well-defined. Another example is in constructing atomic clocks with BEC; information 
about the exact number of atoms is essential in estimating the errors in the frequency. The 
QND schemes proposed are important in that they are all-atom based schemes which utilises 
the inherent collisions directly, rather than using light which may potentially heat or excite 
the ultracold atoms into different quantum states. SC wishes to acknowledge support from 
the Royal Commission for the Exhibition of 1851 and KB the UK EPSRC. 
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FIGURES 
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FIG. 1. Phase shift gained by the displaced condensate after colliding with a stationary condensate at 
the origin. This is given as a ftinction of the nonlinearity, C which is directly proportional to the number of 
atoms in the condensate. The circles represent actual phase shift calculated using a GPE simulation. 
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FIG. 2. Diagram of a possible experimental scheme, (a) Two independent traps (indicated by the 
solid and dash-dot lines) containing condensates in different hyperfine states are prepared and are brought 
towards each other. The probe and reference probe are made by splitting a single condensate, (b) The 
traps are spatially overlapped for the two condensates to interact, (c) One of the traps is then displaced 
and appropriate measurements are made. 
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